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Abstract. We present BeppoSAX LEGS, MEGS, and PDS spectra of eleven X-ray selected BL Lacertae objects. 
Combining these sources with the ones presented elsewhere we have a sample of 21 BL Lacs from the Einstein 
Medium Sensitivity and Einstein Slew Survey. The sample shows strong correlations of several physical parameters 
with the peak frequency of the synchrotron branch of the spectral energy distribution. In particular the peak 
frequency is correlated to the X-ray spectral shape: objects with the peak near to the X-ray band show harder 
and straighter X-ray spectra than those of the low energy peaked sources. This work shows that the recently 
proposed unification scenario for different types of blazars can hold also within the class of high frequency peaked 
BL Lac objects. 
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1. Introduction 

BL Lac objects mark one extreme population of active 
galactic nuclei (AGN). They exhibit intense variability 


(e.g. Wagner & Witzel 1995) and polarization (e.g. Kiihr 


& Schmidt 1990) but do not show strong emission or ab- 


sorption lines in their spectrum (e.g. Stocke et al. |l99l|) . 
In order to understand their physical nature through the 
emission processes involved and their relative contribu¬ 
tion, one has to study their spectral energy distribution 
(SED). 

Determining the SED is a telescope time consuming is¬ 
sue. There exist some well studied bright BL Lac objects - 
like Mkn 421, lES 2344-1-512 and Mkn 501 - for which de¬ 
tailed measurements at several wavelengths are available, 
but to construct a global scenario for BL Lac objects it is 
necessary to also derive the properties of complete sam¬ 
ples. This means investigating fainter objects too. Because 
of its sensitivity and wide energy band, which allows us 
to study simultaneously the BL Lac spectrum from 0.1 to 
100 keV, the BeppoSAX satellite is a powerful tool for this 
aim. 

Here we present a well defined sample of 21 BL Lac 
objects which have been studied using BeppoSAX . In the 
following Section 2 we define the sample; observations and 
data analysis are explained in Section 3. The X-ray data 
are compared to earlier results by the ROSAT satellite in 
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Section 4, and in Section 5 we examine the properties of 
the BL Lac SEDs. The results are discussed in the context 
of unified schemes for BL Lac objects in Section 6. 


2. The sample 

Because strong emission lines are missing in the optical 
spectra of BL Lac objects, this kind of AGN is usually 
selected either in the radio or in the X-ray band. The 
X-ray selected BL Lacs are mainly found in the surveys 
carried out by the EINSTEIN IPC (the EMSS Sample, 
Gioia et al. 1990, and the Slew Survey Sample, Perlman 
et al. |1996D , and by the ROSA T PS PC (ROSA T Al l- 
Sky Survey Sample; Bade et al. 1998 , Beckmann ^000 ). 
The Slew Survey Sample covers the whole high Galactic 
latitude sky, while the EMSS is deeper, but only over 
an area of ~ 800 deg^. By selecting objects with fluxes 
Fjf(0.1 — lOkeV) > 10“^^ergcm“^ sec“^ in the Slew 
Survey and Ex > 4 x 10“^^ergcm“^ sec“^ in the EMSS 
we obtain a sample that combines the advantage of a flux 
limited sample with a wide coverage of the (spectral) pa¬ 
rameter space and thus offers a representative subset of 
the BL Lac population. This allows us to include several 
’’flavors” of BL Lac objects, ranging from the extreme X- 
ray dominated HBL to the borderline objects of interme¬ 
diate BL Lac objects (IBL). The objects we analyse here 
are the second half of the sample, while the first 10 ob- 
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Table 1. The objects of the sample 


3.1. Spectral Analysis 


Name R.A. (2000.0) Dec. Redshift 


lES 

0145-tl38 

01 

48 

29.8 

+ 14 

02 

16 

0.125 

lES 

0323-t022 

03 

26 

13.9 

+02 

25 

15 

0.147 

lES 

0507-040 

05 

09 

38.2 

-04 

00 

46 

0.304 

lES 

0927-t500 

09 

30 

37.6 

+49 

50 

26 

0.186 

lES 

1028-t511 

10 

31 

18.5 

+50 

53 

34 

0.361 

lES 

1118-t424 

11 

20 

48.0 

+42 

12 

10 

0.124 

lES 

1255-t244 

12 

57 

31.9 

+24 

12 

39 

0.141 

lES 

1533-t535 

15 

35 

00.7 

+53 

20 

38 

0.890 

lES 

1544-t820 

15 

40 

15.6 

+81 

55 

04 

7 

lES 

1553-tll3 

15 

55 

43.2 

+11 

11 

20 

0.360 

lES 

1959-t650 

20 

00 

00.0 

+65 

08 

56 

0.047 


jects have been presented in Wolter et al. 
W98. 


(1998), hereafter 


3. BeppoSAX Observations and Data Analysis 


The X-ray observations have been carried out with 
BeppoSAX , a project of the Italian Space Agency (ASI) 
with a participation of the Netherlands Agency for 
Aerospace (NIVR). A detailed description of the entire 
BeppoSAX missi on can be found in Butler & Scarsi ( 1990 ) 
and Boella et al. (|l997a ). The data presented here are from 
the Low Energy Concentrator Spectrometer (LEGS), the 
Medium Energy Concentrator Spectrometer (MECS), and 
the Phoswich Detector System (PDS). 

The LEGS is sensitive in the 0.1-10 keV range with a 
field of view of 37 arcmin diameter (Parmar et al. 1997). 
The LEGS data are useful up to 4 keV only, because the 
response matrix is not well calibrated above this energy 
(Orr et al. 1998). It has an energy resolution of a factor 
~ 2.4 better than the ROSAT PSPC (Brinkmann 1991), 
while the effective area is between ~ 6 and 2 times lower 
at 0.28 and 1.5 keV, respectively. 

The MECS has a field of view of 56 arcmin diameter, 
works in the energy range 1.3-10 keV with an energy res¬ 
olution of ^ 8% and an angular resolution of ^ 0.7 arcmin 
(FWHM) at 6 keV. The effective area at 6 keV is 155 cm^ 
(Boella et ah, 1997b| ). 


The PDS is sensitive in the 13-200 keV band. All pho¬ 
tons within the field of view are counted, with no spa¬ 
tial information (f.o.v. diameter ~ 1.3 deg FWHM). Thus 
other sources in the vicinity of the target can contaminate 
the measurement. 


The 11 objects presented here (for positions and red- 
shifts see Tab. ^ have been observed between May 1997 
and February 1999. The data have been pre-processed at 
the BeppoSAX SDC (Science Data Center) and retrieved 
through the SDC archive. Table || shows the journal of ob¬ 
servations, including exposure times and net count rates 
for the LEGS, MECS, and PDS detectors. 


For observations accomplished in May 1997 all three 
MECS units are available, while after this date only MECS 
units 2 and 3 were in use. The spectra from the MECS 
units have been summed together to increase the S/N. 
None of the sources shows extension neither in the LEGS 
nor in the MECS image. The same reduction process as 
in W98 has been applied to the data, using FTOOLS v4.0 
and XSPEC v.9.0 (Shafer et al. 1991). 

We fitted simultaneously LEGS and MECS 
data, leaving free the LEGS normalization 
with respect to the MECS to account for 
the residual errors in the intensity cross¬ 
calibration. A reliable value for the ratio should 
be LECS/MECS ~ 0.7 ± 0.2 (see the cook¬ 

book for BeppoSAX scientific analysis under 


http://www.asdc.asi.it/bepposax/software/cookbook/ 
and check the NFI flux cross-calibration). For all sources 
where PDS data were available they were also simulta¬ 
neously fitted: only for lES 1255-1-244 there are no PDS 
data and for lES 0145-1-138 the exposure time was too 
short to allow a detection in the PDS. The PDS/MECS 
inter-calibration was handled as a free parameter. In all 
cases the consideration of the PDS part of the spectrum 
did not change significantly the results of the spectral 
analysis, also when fixing the inter-calibration to the 
value of PDS/MECS ~ 0.8 as recommended by the 
BeppoSAX cookbook. In most cases the statistics of the 
PDS was not high, as can be seen from the net counts 
and errors given in Table 

The assumed spectral model is a single-power law plus 
free low energy absorption (arising from cold material with 
solar abundances; Morrison & McCammon 1983). For all 
objects we then checked whether a broken-power law with 
Nh fixed to the Galactic value would give a significantly 
better result (as measured by an F-test value greater than 
95%) than the single-power law fit, and in five cases indeed 
the fit improves. 

We list the best fit parameters for the cases where the 
single power law or broken-power law model gave the bet¬ 
ter fit in Table ^ and Table respectively. Galactic hydro¬ 
gen column values were taken from the Leiden/Dwingeloo 
Survey (Hartmann & Burton 1997) and are listed for refer¬ 
ence. The errors are 90% confidence levels, for three inter¬ 
esting parameters {Nh, cxx, and normalization). Fluxes in 
the 2-10 keV band and also in the 0.5-2.0 keV band (for 
comparison with ROSAT-PSPG fluxes) are given. Also 
listed are the normalization factors of the LEGS rela¬ 
tive to the MECS. For the broken-power law we used the 
LECS/MECS ratio as determined from the single-power 
law fit and fixed Nh to the Galactic value. 

We then also checked if a broken-power law model 
with free fitted absorption would give a better represen¬ 
tation of the spectra. This model increases the number 
of fit parameters and therefore the complexity of the ht. 
Nonetheless we find all free fitted Nh values lying in be¬ 
tween the Galactic value and that from the single-power 
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Table 2. Journal of BeppoSAX observations 


Name 

obs. date 

LEGS 
exp. time [sec] 

LEGS 

net counts 

MEGS 
exp. time [sec] 

MEGS 

net counts 

PDS 

exp. time [sec] 

PDS 

net counts 

lES 0145-bl38 

30-31/12/97 

10576 

73.0 ± 9.5 

12443 

78.7 ± 11.0 

- 

- 

lES 0323-b022 

20/01/98 

6093 

201.6 ± 14.5 

14408 

607.4 ± 27.2 

6845 

256 ± 561 

lES 0507-040 

11-12/02/99 

9116 

441.6 ±21.5 

20689 

1460.2 ± 40.4 

9094 

1515 ±633 

lES 0927-b500 

25/11/98 

8436 

568.4 ±24.3 

22712 

1967.3 ± 46.5 

10129 

692 ± 571 

lES 1028-b511 

1-2/05/97 

4552 

737.9 ±28.1 

12622 

2448.7 ± 50.2 

9484 

2763 ± 718 

lES 1118-b424 

1/5/97 

6027 

236.5 ± 15.6 

9982 

541.3 ± 24.1 

8496 

170 ± 147 

lES 1255-b244 

20/6/98 

2484 

297.1 ± 17.4 

6910 

1037.9 ±33.1 

- 

- 

lES 1533-b535 

13-14/02/99 

8321 

319.4 ± 18.5 

26773 

931.6 ±35.5 

4056 

308 ± 285 

lES 1544-b820 

13/02/99 

8043 

170.7 ± 13.6 

23249 

510.8 ±26.9 

10414 

780 ± 363 

lES 1553-bll3 

5/02/98 

4421 

1179.5 ± 34.5 

10618 

2157.6 ±47.3 

4671 

542 ± 363 

lES 1959-b650 

4-5/05/97 

2252 

423.2 ±20.7 

12389 

3243.4 ± 57.6 

7348 

830 ± 516 


Table 3. BeppoSAX : Best fit results for a single-power law model with free fitted Nh 


Name 

Energy 

Nh 

Nh 

Ftc 


Nm'^ 

xlidof) 

Prob. 


Index ax 

(Gal) 

(Fit) 

2-10 

0.5-2 




0145±138 

1 KA +u.y4 
-L-OU _0.68 

4.59 

10 4 +^^ *' 

_8.9 

0.35 

0.45 

0.86 

0.38 (3) 

77% 

0323±022 

1.^)0 _o.21 

7.27 

31 6 +1^-0 
- 11.5 

2.24 

1.93 

0.81 

0.51 (27) 

98% 

0507-040 

1 14 +012 
-0.11 

7.84 

14 5 

-5.6 

3.88 

2.72 

0.82 

0.73 (69) 

95% 

1028±511 

•1 OQ +0.08 

-0.07 

1.27 

0 '7 +0.8 
-0.7 

10.10 

12.50 

0.66 

0.85 (97) 

85% 

1544±820 

A 1 0 +0.29 
-0.26 

3.70 

20.1 +_8i 

0.95 

2.06 

0.63 

0.68 (30) 

90% 

1959±650 

1 54 +0-08 
-0.08 

0.99 

25.5 til 

12.90 

13.52 

0.65 

0.79 (88) 

92% 


“ hydrogen column density in xlO^°cm“^ 

*’ un-absorbed flux in 10“^^ ergcm“^sec“^ in the 2-10 keV MEGS energy band 
un-absorbed flux in ergcm“^ sec“^ in the 0.5-2.0 keV LEGS energy band 

Normalization of LEGS versus MEGS 


Table 4. BeppoSAX : Best fit results for a broken-power law model applying Galactic Nh values, when significantly 
better representation of the data than the single power-law (according to an F-test). For comparison, also the results 
for single-power law with free fitted Nh are given. The description of the various columns is as in Table ^ 


Name Energy qi Q 2 Eq 

Index ax [keV] 


0927±500 

1.18 

+u.uy 

-0.09 

0.40 

+U.i8 

-0.23 

1.27 

+U.U8 

-0.08 

1.35 

+U.28 

-0.24 

1118±424 

1.57 

+0.16 

-0.16 

1.43 

+0.08 

-0.10 

3.43 

+0.7 

-0.7 

5.11 

+ 1.6 
-2.3 

1255±244 

1.15 

+0.12 

-0.11 

0.61 

+0.19 

-0.37 

1.23 

+0.13 

-0.04 

1.58 

+0.36 

-0.36 

1533±535 

1.57 

+0.15 

-0.14 

0.68 

+0.19 

-0.25 

1.74 

+0.16 

-0.15 

1.40 

+0.30 

-0.27 

1553±113 

1.79 

+0.09 

-0.08 

0.57 

+0.18 

-0.22 

1.85 

+0.09 

-0.08 

1.13 

+0.90 

-0.90 


law and only in one case (lES 1959-1-650) the free Nh 
from the broken-power law is in better agreement with 
the fitted Nh from the single-power law than with the 
Galactic value. Because of that we assumed that in the 
broken-power law model there is no significant low energy 
absorption in the source other than that approximated 
by the Galactic neutral hydrogen column. We thus have 
five free parameters in both models. For the single power 
law with free fitted absorption the spectral slope ax, the 
absorption Nh, the LEGS/MECS and the PDS/MEGS 
normalization, and normalization of MEGS. Although for 
the broken-power-law we added two more free parameters, 
i.e. the break energy Eq and the second spectral index 
02 , we fixed column density Nh to the Galactic value and 
the LECS/MECS normalization to the same value found 
previously for the single power law. All broken-power laws 


N^ 

Nfi 

Ex 

El 

Nm'* 

xlidof) 

Prob. 

(Gal) 

(Fit) 

2-10 

0.5-2 




1.31 

4-3 log 

4.59 

4.59 

0.69 

0.88 (61) 

74% 

2.59 

0 c +1-3 
-1.0 

2.55 

3.65 

0.57 

0.78 (23) 

76% 

1.21 

0 g +1.5 
-1.0 

8.03 

7.81 

0.74 

0.78 (36) 

83% 

1.28 

4 8 

- 1.1 

1.64 

2.89 

0.71 

0.96 (44) 

55% 

3.53 

9.2 

9.37 

19.16 

0.76 

1.05 (89) 

35% 


show a flat slope in the low energy range (oi = 0.4 — 0.7) 
and a steep high energy tail (02 = 1-2 — 1.9) with a break 
energy within the LEGS energy band {Eq = 1.1 —1.6 keV), 
except lES 1118-b424 {Eq = 5.1 keV). 

As mentioned, all sources except lES 1255-1-244 and 
lES 0145-1-138 have been also detected (at > 4 (t) in the 
PDS. In all cases the single- or broken- power law fit give 
a good approximation of the spectrum in the energy range 
0.1 < E < 30keV and even up to 100 keV in the case of 
lES 1028-b511. 

The BeppoSAX data and best fits are plotted in Figure 

0 - 

Before concluding this section, we would like to point 
out the presence of a possible feature in the SAX spec¬ 
trum of lES 0927-1-500. For all models it shows a signifi¬ 
cant excess within the LEGS 1.3 - 1.5 keV energy range. 
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For an added Gaussian line, the fit gives a central energy 
E = 1.4 ± 0.1 keV, corresponding to a rest frame energy 
of ^ 1.7 keV, and a line width of FWHM ~ 200 eV, 
comparable to the energy resolution of the LEGS at this 
energy (F-test probability 99.9%). Since the significance of 
the detection is low and the energy does not have an obvi¬ 
ous physical identification, we refrain from further analysis 
until new data will be available. 


4. Comparison with ROSAT PS PC data 


Four of the 11 sources studied here have also been observed 
during ROSAT PSPC pointed observations and for ten of 
the sources^ data are available from the ROSAT All Sky 
Survey Bright Source Catalogue (RASS-BSC; Voges et al. 


1999). 


In order to compare the ROSAT and 
BeppoSAX results we estimated the best fitting val¬ 
ues for spectral slope and absorption from the ROSAT 
spectra. For the RASS-BSC data we used the hardness 


ratios, a method described by Schartel et al. (|1996D , where 
the hardness ratio is defined as P[R = with H and 
S being the number of counts in the hard and soft energy 
bands; typically two ratios are computed: HRl with 
energy ranges S' = 0.1 — 0.4 keV and H — 0.5 — 2.0 keV, 
and HR2 with S = 0.5 - 0.9keV and H = l.O - 2.0keV 


(Voges et al. 1999). The values for the hardness ratios 
range by definition from -1 for extremely hard to -fl for 
very soft X-ray spectra. The error estimate for the Nh 
and ax values is based on the hardness ratios only, not 
on the photon spectrum itself. Therefore this method 
does not give values, but is able to determine 68%(lcr) 
errors. This is done by exploring the hardness-ratio, 
spectral slope, and Nh parameter space, determining 
the Icr region within it for a given set of parameter 
components. For the pointed observations we instead 
used the standard reduction proce dure as described e.g. 
in Comastri, Molendi & Ghisellini ( 1995 ). 

For all objects we considered two models, namely a 
single-power law with either free fitted absorption or ab¬ 
sorption fixed at the Galactic value. For the pointed ob¬ 
servations the model with free absorption gives acceptable 
results; however the exposure time for the RASS-BSC data 
is quite short and thus we used the values derived from 
the single-power law model with Galactic absorption. The 
results are reported in Tables H and respectively. 


4.1. Comparison between BeppoSAX and ROSAT 
results 

We made a comparison of the fluxes between the ROSAT 
PSPC and the BeppoSAX observations. To do so, we used 
the fluxes in the 0.5-2 keV energy band, that is inside both 
the ROSAT/PSPC and LEGS energy ranges. They vary 

^ The flux of lES 0145-1-138 is too low to be included in 
the RASS-BSC and thus we do not have ROSAT data for this 
object. 



i 1 V 


Fig. 2. X-ray spectral slopes derived from BeppoSAX in 
comparison to those from ROSAT PSPC data. The dashed 
line represents the aseppoSAX = orosat relation. Open 
symbols refer to the RASS-BSC data, while the others 
to ROSAT PS'PG pointed observations. For objects where 
the broken-power law fit to the BeppoSAX spectra gave 
a better statistical result, the lower energy slope ai is 
reported. 


by factors in the range FseppoSAX /-Frosat = 0.3... 2.4 
with a mean ratio FseppoSAX /^rosat = 1-1 ± 0.7. 
Actually, high variability is expected in the X-ray flux o f 
X-ray selected BL Lac objects (e.g. Mujica et al. 1999 ). 
It is interesting that for more radio powerful objects, like 
those belonging to the IJy sample, this ratio is inverted, 
with a measure d me an ratio FReppoSAX /A'rosat = 0.6 
(Padovani et al. 2001 ). Radio selected BL Lac objects have 
flat inverse Compton spectra extending toward high ener¬ 
gies (beyond ~ IkeV), but synchrotron emission, with 
a steep spectrum, may often dominate at soft energies. 
Therefore, the fit of their BeppoSAX LECS/MECS spec¬ 
tra, which sample a much wider energy range toward the 
high energies than the ROSAT PSPC spectra do, may 
be dominated by the flat inverse Compton component. 
Thus the fitted flux at 1 keV is lower than that deter¬ 
mined with fits to ROSAT spectra. This might cause 
the low FReppoSAX /A'rosat ratio computed by Padovani 
et al. This problem does not affect the hard and flat 
soft X-ray spectra of our targets, so that our computed 
FReppoSAX /^rosat ratio is consistent with unity. 

The comparison between the ROSAT and 
BeppoSAX spectral slopes is shown in Figure ^ 
Again we used the spectral indices derived from the 
pointed PSPC observations if available and the RASS- 
BSC spectral slopes otherwise. In fact the best way 
would have been to fit the spectra in the same energy 
ranges for both instruments, as we did for the flux 
values. But we were prevented from doing so by the 
low statistics we would have had. The spectral slope 
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Table 5. ROSAT PSPC: Results for the pointed observations for two models: a single-power law with free fitted and 
with galactic absorption 


Name 

Fx 

NH,Fit 

ElkeV [pJy] 

h 

<^ROSAT 

xlid.ol) 

NH,Gal 

O^ROSAT 

X^(d.o.f.) 

exp.time [sec] 

lES 0323-k022 

4.97 

9.7 ±0.3 

1.57 

1.64 ±0.47 

2.83 (17) 

7.27 

1.39 ±0.25 

7.31 (18) 

25304 

lES 0927-k500 

2.43 

2.6 ±0.4 

0.80 

1.39 ±0.13 

0.90 (17) 

1.31 

0.92 ±0.06 

1.76 (18) 

3831 

lES 1028-k511 

7.78 

1.4 ±0.7 

2.66 

1.51 ±0.03 

3.46 (17) 

1.27 

1.43 ±0.01 

3.65 (18) 

10777 

lES 1118-k424 

1.54 

2.6 ±0.3 

0.48 

1.79 ±0.12 

1.16 (17) 

2.59 

1.79 ±0.03 

1.09 (18) 

6416 

“ un-absorbed flux 

in the 

“hard” PSPC band (0.5- 

-2.0 keV) in 

10 ~^^ erg cm” 

■^sec“’^ 





** spectral index for free fitted Nh 
spectral index for Nh fixed to the Galactic value 


Table 6. ROSAT PSPC: Results for the RASS data for two models: a single-power law with free fitted and with 
galactic absorption 


Name 

Fi 

NH,Fit 

-^IkeV \f^Jy\ 

h 

“ROSAT 

NH,Gal 

Q^rosat 

net counts'^ 

1ES0323±022 

12.3 

18.9 

1.49 

2.58 

+ i.bV 

-1.03 

7.27 

0.94 

301.9 ±24.2 

1ES0507-040 

9.3 

00 

00 

2.54 

0.98 

+0.97 

-0.82 

7.84 

0.82 

274.0 ±23.5 

1ES0927±500 

18.1 

2.2 

2.57 

1.20 

+0.32 

-0.31 

1.31 

0-85 

525.5 ±17.0 

1ES1028±511 

18.0 

1.5 

4.52 

1.39 

+0.23 

-0.22 

1.27 

1 OQ +0-03 
-0.03 

853.9 ±17.7 

1ES1118±424 

4.7 

3.3 

0.80 

1.91 

+0.65 

-0.62 

2.59 

1 04 +008 
-0.08 

139.0 ±8.2 

1ES1255±244 

4.8 

3.8 

0.34 

2.07 

+0.58 

-0.55 

1.21 

1 07 

-0.07 

187.2 ±9.5 

1ES1533±535 

8.0 

2.0 

1.82 

1.23 

+0.23 

-0.23 

1.28 

n QO +0-03 
-0.03 

781.1 ±20.2 

1ES1544±820 

4.0 

5.3 

1.26 

1.77 

+0.51 

-0.50 

3.70 

1 20 +0-09 
-0.09 

252.4 ±13.1 

1ES1553±113 

12.0 

0.7 

3.64 

0.16 

+0.82 

-0.40 

3.53 

1 16 +010 
-0.16 

510.5 ±46.0 

1ES1959±650 

27.9 

15.8 

1.23 

1.76 

+0.44 

-0.21 

0.99 

-0 60 

u.uu _o.oi 

3990.2 ±61.1 


“ un-absorbed flux in the “hard” PSPC band (0.5-2.0 keV) in 10“^^ erg cm“^sec“'^ 
** spectral index for free fitted Nh 
spectral index for Nh fixed to the Galactic value 
PSPC 0.5-2.0 keV band 


values vary between the different observations, but the 
average difference is anosAT - ctBeppoSAX = 0.0 ± 0.4 
(if we consider only the objects with PSPC pointed 
observations we get anoSAT - ctBeppoSAX = 
0.17 ± 0.07, while for the RASS-BSC data 
o-rosat — oiBeppoSAx = —0.1 ± 0.5). While these 

spectral indices are consistent with each other it is 
remarkable that the broken-power law model applied 
to the BeppoSAX data shows in the low energy band 
{E < 1.5 keV) a slope flatter than that of the ROSAT 
PSPC (0.1 - 2.4 keV). The difference might be a hint 
of more complex X-ray spectra or could result from the 
ROSAT PiSPC calibration problems discussed by Iwasawa 
et al. ( 1999D and Barcons et al. (200C), who showed the 
X-ray spectral slopes measured by the ROSAT PSPC 
to be significantly steeper than those from other X-ray 
missions by Aax ~ 0.4. PSPC spectra steeper than the 
BeppoSAX ones have also been found for the IJy sample 
(Padovani et al. 2001). 


5. Spectral Energy Distribution 

To determine the broadband SED of the 11 objects we 
used radio data at 1.4 GHz from the VLA surveys NVS S 
(Condon et al. 1998 ) or FIRST (White et al. 1997 ). 
Optical data were taken from the literature and for some 


objects dete rmined using the Calar Alto 1.23m telescope 
(Beckmann ^OOOa ). A campaign for simultaneous optical 
da ta has also been performed for a few objects (Villata et 
al. 20001) . Variability in the optical band is not expected to 
be large: all objects presented here are X-ray dominated 
objects (i.e. aox < 1-2, see Table |^) w hich show only 
small optical variability (e.g. Villata et al. 2000, Mujica et 
al. 1999, Januzzi et al. 1994). With these data available 
we also computed the broad band spectral indices (be¬ 
tween 1.4GHz, 4400 A and 1 keV), which are reported in 
Table |. 


The following step consisted in parameterizing the de¬ 
rived SED. 


Several authors pointed out that the synchrotron 
branch of the BL Lac SED can be well approximated 
by a parabolic spectral shape (cf. Landau et al. 


1986 


Comastri et al. 1995, Sambruna et al. 1996, Fossati et al. 
[1998 , W98). We thus applied a parabolic fit to the SED 
in the log v — log plane - parameterized as log = 
a ■ (log h')'^ + b- log v + c- and determined the peak position 
(t^peak) and the total flux/luminosity of the synchrotron 
component. The resulting Vpeak and two-point indexes are 
listed in Table ^ We also computed these quantities for 
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the objects reported in W98 q In the following discussion 
the whole sample of 21 BL Lac objects will be considered. 

While the parabolic fits just described give a useful 
general parameterization of the SED, we have also con¬ 
sidered all the available data and reproduced the SED of 
all sources of our sample with model dependent represen¬ 
tations. 


5.1. The blazar spectral sequence 


Fossati et al. ( 1998D proposed a simple phenomenologi¬ 
cal description of the average SED of blazars, based on 
the their bolometric observed luminosity which is in turn 
assumed to be well traced by the radio one. The latter 
in fact was found to correlate (and thus define) both the 
peak frequency of the synchrotron spectrum and the rel¬ 
ative importance of the inverse Compt on vs synchrotron 
powers. More recently, Donato et al. ( 2001 ) have revis - 
ited the original parameterization by Fossati et al. ( 1998|) , 
but (only) for objects below a radio luminosity of 10 "^^ erg 
s“^, assuming a slightly different relation between the ra¬ 
dio power, the synchrotron peak frequency and the i nverse 
Compton luminosity. Thus in the Donato et al. ( ^001 ) 
model low power sources are assumed to have comparable 
luminosities in the synchrotron and self-Compton spectral 
components. We have then applied this parameterization 
to the objects in our sample, and the resulting ‘fits’ are 
shown as dashed lines in Fig. ^ 


5.2. Homogeneous Synchrotron Self-Compton models 

Finally we have reproduced the energy distributions 
adopting a more ‘physical’ approach by using an homo¬ 
geneous (one-zone) synchrotron self-Compton model (this 
is a “one-zone version” of the model described in detail 


in Spada et al. 2001): it assumes that the source is cylin¬ 


drical, of size R and width AR' = R/T (in the comoving 
frame, here T is the bulk Lorentz factor). The model is 
aimed at reproducing the spectrum originating in a lim¬ 
ited part of the jet, the one thought to be responsible for 
most of the emission. This region is necessarily compact, 
since it must account for the fast variability shown by all 
blazars, especially at high frequencies. Therefore the ra¬ 
dio emission from these compact sources is strongly self- 
absorbed, and the model cannot account for the observed 
radio flux. 

Sources in our sample are all BL Lacs of relatively 
low power with no signs of strong broad emission lines 


They slightly differ from the values presented there due to 
the different reference energies: our aox, ctro and On are lower 
by ~ 0.17, ~ 0.05, and ~ 0.09 than those presented in W98. 
Furthermore after the publication of W98 two more redshifts 
for objects of the sample have been determined, lES 0502-1-675 
{z = 0.314, Scar pa et al. |l999| ) and lES 1517-b656 [z = 0.702, 
Beckmann et al. 1999). Thus we lack only redshift information 
for lES 1544-1-820, for which we assume 2 = 0.2. 

® Here not intended in a statistical sense, but only as approx¬ 
imate representations of the SEDs 


in their spectra. Consequently, we have neglected the role 
of external seed photons for the inverse Compton (IC) 
process. 


Relativistic electrons between a minimum and a maxi¬ 
mum energy are injected throughout the source for a lim¬ 
ited amount of time tinj (mimicking a flare), which we 
set equal to tinj = AR'/c. High energy electrons whose 
radiative loss timescales are shorter than tjnj reach a sta¬ 
tionary distribution, while lower energy electrons retain 
the original spectrum. The electron energies for which the 
cooling timescale is equal to tinj are denoted by 7 c. The 
particle distribution N{y) is assumed to have the slope 
n [iV( 7 ) oc 7 “"] above 7 c, while below this value there 
can be two cases. We in fact assume that the particle dis¬ 
tribution derives from a continuous injection of particles 
between 7 min and 7 niax, and the slope below 7 c depends on 
whether 7 c is greater or smaller than 7 niin- If 7 c > 7 min, we 
have N{y) oc 7 “"+^ between 7 inin and 7 c and N{y) oc 7 “^ 
below 7 niin- Alternatively, if 7 c < 7 min, then N{y) oc 7 “^ 
between 7 c and 7 min- We further assume that, below the 
minimum between 7 c and 7 min, N{y) oc 7 “^. According 
to these assumptions, the random Lorentz factor 7 peak of 
the electrons emitting most of the radiation (i.e. emit¬ 
ting at the peaks of the SEDs) becomes equal to 7 c and 
is determined by the importance of radiative losses and 
can assume values within the range 7 min“ 7 max- Its value 
is listed in the last column of Table The source is as¬ 
sumed to emit an intrinsic luminosity L' and is assumed 
to be observed with the viewiM angle 9. All these input 
parameters are listed in Table ||. 


As can be noted from it, the intrinsic luminosities, the 
source dimensions, the bulk Lorentz factors and viewing 
angles, and the magnetic field values are very similar for 
all sources, in agreement with the fact that all objects 
belong to the same “flavor” of blazars. The resulting fits 
are shown in Fig. ^ as solid lines. 


Let us compare these representations of the SEDs and 
in particular the resu lts of the parameterization according 
to the Fossati et al ( 1998 )’s scenario with those inferred 
from the homogeneous SSC model. It turns out that the 
former tends to overestimate the Compton 7 -ray emis¬ 
sion (by construction its power is never less than the syn¬ 
chrotron one), but it agrees with the existing upper limits 
in the GeV and TeV bands in all cases but lES 1553-1-113 
and lES 1959-1-650. We stress that both these parameter- 
izations are not meant to describe accurately the SED of 
specific sources, but only the average SED of sources of 
equal total power. Keeping this in mind, we can consider 
the difference between the two theoretical spectra for each 
source as a measure of the uncertainty associated with the 
‘theoretical’ description of the SED. 


5.3. Derived parameters and correlations 

We can now consider possible trends in the properties 
of the SED and examine their consistency with the pre- 
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Table 8. Input parameters of the homogeneous synchrotron self-Compton model 


Name 

erg 

R 

cm 

B 

G 

r 

e 

n 

'ymin 

Tpeafe 

0145-kl46 

6 .0x10'*“ 

8 .0x10*° 

1.0 

11 

5.0 

3.8 

600 

27000 

0323-k022 

1 .0x10*^ 

1 .0x10*° 

0.9 

11 

3.7 

3.5 

800 

27000 

0507-040 

2.5x10*^ 

1 .0x10*° 

1.2 

15 

3.5 

3.3 

1200 

20000 

0927-k500 

6 .0x10'*^ 

1.5x10*° 

0.7 

14 

4.8 

3.4 

400 

30000 

1028-k511 

2.7x10'*^ 

1 .0x10*° 

1.0 

14 

3.0 

3.2 

800 

26000 

1118-k424 

2 .0x10'*^ 

1 .0x10*° 

1.5 

15 

5.0 

3.3 

300 

13000 

1255-k244 

6.5x10*^ 

1 .0x10*° 

1.0 

12 

5.0 

3.3 

500 

18000 

1533-k534 

1.5x10*^ 

2 .0x10*° 

1.0 

15 

2.6 

3.2 

200 

12000 

1553-kll3 

2.5x10*^ 

3.0x10*° 

0.7 

15 

2.5 

3.55 

300 

13000 

1959-k650 

8 .0x10'*° 

1 .0x10*° 

1.2 

13 

4.0 

3.6 

500 

19000 


Table 7. Derived quantities: Two-point broad band spec¬ 
tral indices“, X-ray luminosities and peak frequencies (for 
the 11 objects from this paper and for the 10 BL Lacs 
from W98). 


Name 

Qiox 

Qiro 

CTrx 

logL^f 

log(i^peafc) 

lES 0145-kl38 

1.19 

0.41 

0.65 

43.42 

14.49 

lES 0323-k022 

1.07 

0.34 

0.57 

44.38 

14.89 

lES 0507-040 

0.69 

0.52 

0.58 

45.26 

18.05 

lES 0927-k500 

0.81 

0.36 

0.50 

44.89 

16.34 

lES 1028-k511 

0.81 

0.33 

0.48 

45.86 

16.08 

lES 1118-k424 

0.91 

0.32 

0.51 

44.30 

15.49 

lES 1255-k244 

1.14 

0.15 

0.46 

44.88 

15.09 

lES 1533-k535 

0.77 

0.37 

0.51 

46.05 

15.77 

lES 1544-b820° 

0.71 

0.47 

0.55 

44.11 

16.15 

lES 1553-kll3 

0.81 

0.43 

0.56 

45.89 

15.58 

lES 1959-k650 

1.05 

0.31 

0.54 

44.11 

15.00 

MS 0158.5-k0019 

0.72 

0.39 

0.50 

45.12 

16.95 

MS 0317.0-kl834 

0.64 

0.42 

0.49 

45.11 

18.51 

lES 0347-121 

0.67 

0.34 

0.44 

45.05 

16.97 

lES 0414-b009 

0.76 

0.42 

0.53 

45.59 

16.25 

lES 0502-b675 

0.59 

0.34 

0.42 

46.00 

18.10 

MS 0737.9-k7441 

0.84 

0.39 

0.54 

44.93 

15.72 

lES 1101-232 

0.72 

0.39 

0.49 

45.83 

17.37 

lES 1133-k704 

1.04 

0.43 

0.62 

43.69 

14.96 

MS 1312.1-4221 

1.04 

0.25 

0.50 

44.85 

15.15 

lES 1517-k656 

0.75 

0.32 

0.47 

46.55 

16.17 

“ derived from source fluxes at 1 keV, 4400 A, and 1.4 GHz, 


resp. 

*’ luminosity in the 2-10 keV band (MEGS) in [erg/sec] 
assuming a redshift z = 0.2 


dictions of the blazar unifying scenario discussed above 
(Fossati et al 1998, Ghisellini et al 1998| ). 


The first point we want to stress is that we find a 
strong correlation between aox and the value of the peak 
frequency. The correlation is well represented by two lin¬ 
ear regressions with a break around 3 x 10^® Hz as shown 
in Figure (Note that small values of Uox refer to X-ray 
dominated objects while more optical dominated objects 
would have aox > !■) Indeed this result is similar to that 



Fig. 4. X-ray dominance (aox) versus the peak frequency 
of the synchrotron branch. The more X-ray dominated 
objects have a higher peak frequency. The lines show two 
linear regression fits with a break at logVpeak = 16.5. 


found by Fossati et al. ( 1998 ) when studying the depen¬ 
dence of Uro and Urx on the peak frequency. 

A second interesting trend appears to be present be¬ 
tween the intrinsic absorption ANh = Nh^fu — Nn^Cai 
and the X-ray spectral slope ax- For all fits the value of 
the absorbing column was larger than the Galactic value. 
Even though error bars are large and part of the cor¬ 
relation might be spurious, there is indication of higher 
ANh for steeper X-ray spectra (Figure ||). If this trend 
were due to intrinsic absorption we should also find a 
correlation with the X-ray luminosity, which instead is 
not observed. Alternatively, in a single power law model 
the absorption can mimic an intrinsic curvature. Indeed 
a curvature of the X-ray spectra can easily account for 
the ANh vs ax behavior: near the synchrotron peak the 
spectra are harder and straighter than in the steep de¬ 
cay at frequencies higher than the peak. And in fact, the 
broken-power law model shows Nh values in agreement 
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I 


log(NH.n, 

Fig. 5. Logarithmic “intrinsic” absorption {Nh^fu — 
Nn.Gai in 10^°cm“^) versus X-ray spectral slope. There is 
a trend for steeper spectra to correspond to higher “intrin¬ 
sic” absorption. Filled symbols represent the cases where 
the single-power law fit was a better approximation to the 
BeppoSAX spectra than the broken-power law. 



Fig. 6. X-ray spectral slope versus peak frequency. 


with the Galactic hydrogen column density (when the 
BeppoSAX spectrum has enough counts to give a good 
statistics for this fit), and in all cases the value for a free 
fitted Nh for the pointed observations is in between the 
Galactic and the free fitted value from the single-power 
law. Furthermore we also find that a higher peak frequency 
is associated with a flatter spectral slope (Fig. ||), a trend 
that again can be accounted for by the above behavior: 
when the peak frequency is rising the X-ray band is lo¬ 
cated near the maximum of the synchrotron branch and 
this gives rise to flatter spectra. 


A further correlation appears to be present between 
the X-ray dominance and the X-ray luminosity (Figure 
Using a linear regression we find that it can be described 
as log Lx = 47.7 —3.2aox with a correlation coefficient of 
0.7. This cannot be due to selection effects (e.g. missing 
optical faint counterparts of the weaker X-ray sources). 
In this case there should be also a correlation between 
the X-ray fluxes and the X-ray dominance, but this is not 
detectable. This trend was already reported by Beckmann 
( 1999| , [2000 ) using a complete sample of BL Lac objects 
based on the RASS-BSC. 

The same effect of X-ray dominated objects appearing 
to have X-ray luminosities higher than the intermediate 
BL Lacs should also be seen by comparing the peak fre¬ 
quency of the synchrotron branch with the X-ray lumi¬ 
nosity. Within the scenario for the spectral properties of 
blazars discussed so far it should also be expected that 
as the objects considered here have similar radio power 
(i.e. total luminosity), their X-ray luminosity is bound 
to be mostly affected by the location of the synchrotron 
peak, higher energies corresponding to higher luminosi¬ 
ties. However, contrary to these predictions, there does 
not seem to be a correlation between the peak frequency 
and the total luminosity of the synchrotron branch (as 
determined through the parabolic fit to the radio, optical 
and up to the X-ray data as described above). This might 
indicate that when objects of the same class are considered 
the range in total power is too small to show such a trend. 
The range of luminosities which is covered by the sample 
presented here is indeed small: the radio luminosity ranges 
from Lr = 1.2 x lO^^W/Hz to Lr = 1.6 x lO^^W/Hz, 
while the spread in the X-ray luminosity is one decade 
larger (1.4 x 10^®W/Hz < Lx < 1.7 x lO^^W/Hz). For the 
radio, optical, and X-ray energy band studied here, we find 
an average oflogvLR = 41.1 ±0.6, logi^Lopt = 44.7±0.6, 
and log FLx = 45.0 ± 0.8, which are consistent with t he 
average SEDs of HBL as shown by Fossati et al. (1998). 


6. Discussion and Conclusions 

The (physical) parameters of the objects studied here seem 
to partly correlate with the value of the peak frequency 
of the synchrotron branch of the SED. The X-ray se¬ 
lected sample considered in this work comprises objects 
whose synchrotron peaks are located near/in the X-ray 
band (high frequency peaked BL Lac, HBL; Padovani & 
Giommi 1996). Those of radio selected BL Lac objects 
tend instead to be located in the optical or infrared bands 
(low frequency peaked BL Lac, LBL), as is the case of 
sources belonging to the 1-Jy BL Lac sample (Stickel et 
al. 1991) whose peak frequencies can be as low as ~ lO^^Hz 
(Padovani et al. 2001). 

We find a correlation between the peak frequency and 
the spectral slope as well as between the spectral slope 
and the flattening of the soft X-ray spectrum (parameter¬ 
ized as “intrinsic” absorption). These trends are consistent 
with the properties of the SED of HBLs: when the peak 
of the SED moves into the X-ray region the X-ray spectra 
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Fig. 7. X-ray luminosity for the BeppoSAX MEGS band 
(2-10 keV) versus X-ray dominance {aox)- Only for lES 
0145-1-138 (bottom right) the error on the luminosity 
based on the count rate statistics is larger than the plotted 
symbol. 


energy part of the IC component (Padovani et al. ^001 


W98). Interestingly, no flattening associated with the lat¬ 
ter component is observed up to several tens of keV (as 
the PDS data show that the X-ray power-law can be ex¬ 
tended up to lOOkeV), implying that the IC seems to 
start emerging at frequencies v > 10^®Hz. This is fully 
consistent with the assumed SED for HBL as inferred from 
the sequence scenario (Fossati et al. ( 1998| ). 

In summary the BeppoSAX spectral survey shows that 
the X-ray properties of X-ray selected BL Lac objects 
are in good global agreement with the unified model 
for blazar, which ascribes the differences among blazars 
mostly to the location of the peak frequencies of the syn¬ 
chrotron and inverse Compton spectral components (e.g. 
Padovani & Giommi 1996, Ghisellini et al. 1998| ). 

Even though the correlation of peak freq uency with to¬ 
tal luminosity as described by Fossati et al. (|199^ ) cannot 
be found in this sample, we find that the spectral prop¬ 
erties are fully consistent with objects having a spread in 
the position of the peak of the synchrotron branch around 
the X-ray band. The sample as a whole represents only 
a part of the blazar population. Comparing the luminosi¬ 
ties in the different energy bands studied here, the objects 
match the average SEDs of HBL as shown by Fossati et 


are expected to become harder and straighter as the peak 
frequency rises. 

This behavior also accounts for the relative high X- 
ray luminosities observed in X-ray dominated BL Lacs 
(although LBL have higher total bolometric luminosities 
than HBL their X-ray emission corresponds to the energy 
band between the synchrotron and IC components, e.g. 
Fossati et al. 1998 ). Although within our sample it has 
not been possible to determine that HBL have lower bolo¬ 
metric luminosity (as a large fraction of the bolometric 
luminosity is expected to be emitted in the IC branch, 
which extends to energies outside the frequency range we 
covered), we can assume - as already mentioned - that 
the radio luminosi ty is a good tra cer of the bolometric 
one (Fossati et al. 1998, Ghisellini 1999). The lack of a 
visible trend between such luminosities and the peak fre¬ 
quencies has then to be ascribed to the small luminosity 
range spanned by the HBL sample in comparison with 
the intrinsic dispersion of the peak-luminosity relation (as 
expected e.g. for a distribution of the observing viewing 
angle which tends to produce an opposite peak-luminosity 
trend). 

If we consider this, the absence of correlation between 
the total luminosity in the synchrotron branch and the 
peak frequency will then result in higher bolometric lumi¬ 
nosities for the LBL in comparison to the HBL. This is 
consistent with the fact that LBL show strong IG com¬ 
ponents compared to the synchrotron branch, while HBL 
emit similar luminosities in the IC and the synchrotron 
branch. 

As the sample discussed here only includes HBL ob¬ 
jects, we would expect the X-rays to be dominated by the 
higher energy part of the synchrotron and/or the lower 


The energy coverage provided by BeppoSAX toward 
the high energies did not allow us to directly examine 
the link between the position of the peaks and the bolo¬ 
metric emitted power, as postulated by such model. Such 
study will be however possible with observations with sen¬ 
sitive instruments and good spectroscopy resolution above 
~ 300 keV as provided by e.g. the INTEGRAL mission, 
which will cover the energy range of the spectral change 
from the synchrotron to the IC dominated emission and 
possibly of the peak of the IC component. 
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Fig. 3. Spectral energy distribution of our sources. The solid lines correspond to the one-zone homogeneous syn¬ 
chrotron and inverse Compton model calculated as explained in the text, with the parameters listed in Table The 
dashe d lines represent instead the spectrum predic ted b y the phenomenological SED description by Fossati et al. 
( |1998[ ) (with the changes proposed by Donato et al. 2001 ). The iJeppoSAX data and simultaneous optical data of our 
observing campaigns are indicated by the filled bow-ties and symbols. 
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Fig. 3. (continued) 






